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We propose a simplified version of the Kitaev's surface code in which error correction requires 
only three-qubit parity measurements for Pauli operators XXX and ZZZ. The new code belongs to 
the class of subsystem stabilizer codes. It inherits many favorable properties of the standard surface 
code such as encoding of multiple logical qubits on a planar lattice with punctured holes, efficient 
decoding by either minimum-weight matching or renormalization group methods, and high error 
threshold. The new subsystem surface code (SSC) gives rise to an exactly solvable Hamiltonian with 
3-qubit interactions, topologically ordered ground state, and a constant energy gap. We construct 
a local unitary transformation mapping the SSC Hamiltonian to the one of the ordinary surface 
code thus showing that the two Hamiltonians belong to the same topological class. We describe 
error correction protocols for the SSC and determine its error thresholds under several natural 
error models. In particular, we show that the SSC has error threshold approximately 0.6% for the 
standard circuit-based error model studied in the literature. We also consider a model in which 
three-qubit parity operators can be measured directly. We show that the SSC has error threshold 
approximately 0.97% in this setting. 



I. INTRODUCTION 

Quantum error correcting codes are vital ingredients in 
all scalable quantum computing architectures proposed 
so far. By actively monitoring and correcting errors, the 
encoded quantum states can be protected from noise up 
to any desired precision provided that the error rate of 
elementary quantum operations is below certain constant 
value known as the error threshold. 

Topological codes such as the surface code family [IH3] 
have received considerable attention lately due to their 
several attractive features. First, the quantum hardware 
envisioned in the surface code architecture consists of a 
2D array of qubits with controlled nearest-neighbor in- 
teractions and a local readout. In principle, it can be 
implemented using the Josephson junction qubits tech- 
nology [3] . Surface codes feature an error threshold of at 
least 1% [5] which is one of the highest thresholds among 
all studied codes. Secondly, encoded Clifford gates such 
as the CNOT gate can be implemented efficiently by 
the code deformation method [SHE] which requires only 
a mild overhead in space and time. The error rate of 
encoded gates decreases exponentially with the lattice 
size [9]. Thirdly, the surface codes can be decoded effi- 
ciently using Edmonds's minimum weight matching algo- 
rithm [TU] or renormalisation group methods [TTHT^] . 

Although the surface code is among the best code can- 
didates, a promising direction for improvements has been 
recently identified by Bombin [14] who proposed topo- 
logical subsystem codes [15 . A subsystem code [151 HZ] 
can be viewed as a regular stabilizer code in which one 
or several logical qubits do not encode any information. 
The presence of unused logical qubits, known as gauge 
qubits, simplifies eigenvalue measurements of multi-qubit 
stabilizers — such as the plaquette and star operators of 
the surface code — which are required for error correction. 
Consider as an example the simplest 4-qubit code with 



two stabilizers S x = X m and S z = Z®\ It encodes 
two qubits with logical Pauli operators X l = X1X2, 
Z L = Z X Z Z and X G = XtX 3 , Z G = Z 1 Z 2 . If only the 
first logical qubit is used to encode information, the syn- 
drome (eigenvalue) of S can be determined indirectly 
by measuring eigenvalues of the unused logical opera- 
tors X G and S X X G = X 2 X 4 . Multiplying the mea- 
sured eigenvalues together yields the desired eigenvalue 
of S x . The syndrome of S z is determined similarly by 
measuring eigenvalues of the unused logical operators Z G 
and S z Z G = Z3Z4 followed by multiplication of the out- 
comes. Hence the full syndrome extraction requires only 
two-qubit parity measurements and simple classical post- 
processing. The unused logical operators that need to be 
measured in order to extract the syndrome of all stabi- 
lizers are usually referred to as gauge generators, see [17] 
for the general theory of subsystem codes. 

A simplified syndrome readout offered by subsystem 
codes has its own costs. In any practical settings, eigen- 
values of individual gauge generators can only be mea- 
sured with a finite accuracy. As one multiplies together 
measured eigenvalues, errors tend to accumulate render- 
ing the inferred syndrome bit unreliable. This strongly 
limits the class of candidate subsystem codes for the 
topological fault-tolerant architecture. First, a suitable 
code must have local gauge generators, ideally, 2- or 3- 
qubit Pauli operators acting on nearest-neighbor qubits. 
This ensures that the syndrome readout requires only 
local measurements. To avoid accumulation of measure- 
ment errors, a suitable code must also have low-weight 
stabilizers. More precisely, each stabilizer must be com- 
posed of only a few gauge generators. The latter require- 
ment leaves out many interesting families of codes, such 
as the 2D Bacon-Shor codes [TB] and random 2D sub- 
system codes discovered in [18]. In contrast, subsystem 
color codes found in [14] have 2-qubit gauge generators 
while stabilizers act on either 6 or 18 qubits. These codes 
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were shown to have a constant error threshold of at least 
2% under depolarizing noise assuming noiseless syndrome 
readout 19, 20 . Unfortunately, subsystem color codes do 
not inherit favorable properties of the surface code such 
as encoding of multiple logical qubits |14j on a planar 
lattice required for the code deformation method. 

In the present paper we propose a subsystem version of 
the standard surface code on the regular square lattice. 
Each plaquctte of the lattice carries one gauge qubit and 
a pair of weight-6 stabilizers of type X® 6 and Z® 6 , see 
Fig. [T] The code has 3-qubit gauge generators of type 
XXX and ZZZ which makes it suitable for architec- 
tures where direct 3-qubit parity measurements in the 
X- and Z-basis are available. A promising proposal for 
implementing 3-qubit parity measurements in Josephson 
junction qubits has been recently made by DiVincenzo 
and Solgun [21]. By analogy with the surface code, the 
new subsystem surface code (SSC) can encode multiple 
logical qubits on a planar lattice with punctured holes. 
We describe error correction protocols for the SSC and 
determine its error thresholds under several natural er- 
ror models. First, we study the so called "code capacity" 
model where each qubit is subject to independent bit-flip 
and phase-flip errors with rate p, while syndrome mea- 
surements are noiseless. By relating the optimal error 
correction to the phase transition in the random-bond 
Ising model on the honeycomb lattice, we show that the 
threshold error rate is p ~ 7%. Secondly, we study the 
circuit-based error model where the syndrome readout is 
simulated by noisy quantum gates, single-qubit measure- 
ments, and ancilla preparations. Each operation can fail 
with a probability p, see Section VI for details. Monte 
Carlo simulation suggests that the SSC has error thresh- 
old p c « 0.6% for the circuit-based error model. Finally, 
we consider a model in which 3-qubit parity can be mea- 
sured directly. We show that the SSC has error threshold 
approximately 0.97% for this direct parity measurement 
model. 

The new code also gives rise to an exactly solvable 
Hamiltonian with 3-qubit interactions which has a topo- 
logically ordered ground state and whose excitations are 
non-interacting abelian anyons. The exact solvability of 
the model stems from a peculiar commutativity struc- 
ture of the model's Hamiltonian. We show that the set 
of all relevant 3-qubit interactions can be partitioned into 
small clusters such that interactions from different clus- 
ters pairwise commute. In contrast, it was shown re- 
cently by Aharonov and Eldar [22] that topological order 
cannot be realized by 3-local Hamiltonians in which all 
interactions pairwise commute. We also construct a local 
unitary transformation U that maps the 6-qubit stabiliz- 
ers of the SSC to the plaquette and star operators of the 
Kitaev's surface code on the square lattice. The gauge 
generators XXX and ZZZ are mapped to single-qubit 
X and Z operators. Loosely speaking, the map U decou- 
ples gauge qubits from the surface code qubits. 

The rest of the paper is organized as follows. Sec- 
tions II|III introduce a subsystem version of the toric 



code with 3-qubit gauge generators, the corresponding 
exactly solvable Hamiltonian and discuss its connection 
to the ordinary toric code. Extension to the planar geom- 
etry is given in Section [TV] that defines subsystem surface 
codes. The concept of a virtual lattice which is crucial 
for understanding our error correction protocols is intro- 
duced in Section [V] This section also discusses error cor- 
rection for the idealized settings when syndrome readout 
is noiseless. Error correction protocols for the circuit- 
based syndrome readout and numerical simulations are 
presented in Section |VI| Finally, Section |VII| focuses on 
a model in which direct 3-qubit parity measurement are 
available. 



II. A SUBSYSTEM TORIC CODE 

We begin by introducing a subsystem version of the 
toric code that encodes two logical qubits 0Q. Extension 
to planar lattices with a boundary will be described in the 
next section. The code is defined on the regular square 
lattice of size Lx L with periodic boundary conditions. It 
contains L 2 vertices, 2L 2 edges, and L 2 plaquettes. We 
place a qubit at every vertex and at the center of every 
edge of the lattice. Hence there are in total n = 3L 2 
code qubits. For each plaquette p define weight-6 Pauli 
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FIG. 1: Subsystem toric code. Qubits live at vertices and 
centers of edges of the regular square lattice. Opposite sides 
of the lattice are identified. Left: Four types of triangles and 
the corresponding triangle operators G(T). Triangle oper- 
ators that belong to different plaquettes pairwise commute. 
Right: Stabilizer operators Sp (top) and Sp (bottom). Sta- 
bilizers are analogues of the plaquette and star operators of 
the standard toric code. Triangle operators commute with 
stabilizers. Eigenvalue of any stabilizer can be determined by 
measuring eigenvalues of individual triangle operators. For a 
lattice of size Lx L the code has parameters [[3L 2 , 2, L]]. This 
should be compared with the standard surface code which has 
parameters [[2L 2 , 2, L]]. 
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operators S* and Sp as shown on Fig. [T] (right). One 
can easily check that Sp and Sg commute with each 
other for all p and q. Let S be the abelian group gener- 
ated by 2L 2 operators Sp and Sp . It defines a quantum 
stabilizer code with a codespace C spanned by n-qubit 
states ip invariant under S, that is, if) G C iff Sp ■ ip — if) 
and Sp ■ if) = if) for all p. A simple algebra shows that 
Y\ p Sp = I and Yip Sp = I, where the product is taken 
over all plaquettes of the lattice. Furthermore, since each 
qubit belongs to exactly two stabilizers Sp and two sta- 
bilizers Sp , these are the only dependencies among the 
generators of S. This shows that S has s = 2(L 2 — I) 
independent generators. The standard stabilizer formal- 
ism [23] then implies that S is a stabilizer code encoding 
k'=n-s = L 2 + 2 qubits, that is, dim(C) = 2 k ' . 

We shall now divide the k' encoded qubits into g = L 2 
gauge qubits and k = 2 logical qubits. Let u be any 
vertex of the lattice and /, g be a pair of orthogonal edges 
incident to u. The triple (u, /, g) will be referred to as 
a triangle. Note that the lattice has four non-equivalent 
types of triangles, see Fig. [I] We shall say that a triangle 
T = (u, /, g) is north-west (NW) if u is at the north-west 
corner of the plaquette formed by / and g. Similarly one 
defines north-east (NE), south-west (SW), and south-east 
(SE) triangles. Define triangle operators 



G{T) 



X u X f X g if T is SW or NE triangle, 
Z u ZfZ g if T is SE or NW triangle, 



(1) 



see Fig. [T] Here subscripts indicate qubits acted upon by 
the Pauli operators X and Z. Note that triangle oper- 
ators that belong to different plaquettes commute with 
each other. 

Any stabilizer can be expressed as a product of two 
triangle operators using identities 



dp — U(l p )U{l p ), 
oZ _ r , i r rSE\ri(^rNW\ 



(2) 



see Fig.fi] Here T p w , T* E , T p SM/ , and T^ E are triangles 
-pe N 



p ' p 

of type TN W, NE, SW, SE respectively that belong to a 
plaquette p. 

Let us now show that triangle operators commute with 
all stabilizers, thus being logical operators of the code S. 
Consider any stabilizer, say, Sp and any triangle oper- 
ator G{T) of Z-type. If T does not belong to the pla- 
quette p then Sp commutes with G(T) because triangle 
operators from different plaquettes always commute. If 
T belongs to p then G(T) anti-commutes with both X- 
type triangles forming S* , that is, G(T) commutes with 
Sp . A similar argument shows that Z-type stabilizers 
commute with X-type triangle operators. 

The above observations show that we can choose g = 
L 2 pairs of logical operators for the code S as 



Xp^GiT^) and Z p = G{T; 



SE\ 



(3) 



where p runs over all plaquettes of the lattice. We shall 
treat encoded qubits defined by X p and Z p as gauge 



qubits encoding no useful information because the corre- 
sponding logical operators have very small weight. Hence 
each plaquette of the lattice carries one gauge qubit. As 
we will show in Sec. |III| it is possible to completely dis- 
entangle these gauge qubits from the code with a depth 4 
local quantum circuit, leaving behind the usual toric code 
on 2L 2 qubits and L 2 ancillary qubits that are decoupled 
from the code. 

Recall that the code S has k' — L 2 + 2 encoded qubits. 
This leaves k = k' — g = 2 logical qubits which have not 
been identified yet. Let T and A be the set of all qubits 
lying on some fixed horizontal and some fixed vertical 
line of the lattice respectively, see Fig. [2} Note that |T| = 
|A| = 2L. Define 

X 1 = ]JX j , Z 1 = ]JZ j (4) 

jer jga 



and 



X 2 =HXj, Z 2 = ]JZ j . (5) 




FIG. 2: Non-contractible loops on the toric code lattice. Each 
loop contains 2L qubits. 

Any triangle operator commutes with Xj and Zj be- 
cause triangles share or 2 qubits with V and A. In 
addition, since T and A are even sets of qubits overlap- 
ping on exactly one qubit, one has commutation rules 
XiZj = (— l) Si ' j ZjXi. Hence one can view Xi,Zi and 
X 2 , Z 2 as logical X and Z operators on the two remain- 
ing logical qubits of the code S. 

The minimum distance d of a subsystem code is defined 
as the minimum weight of a Pauli error E that commutes 
with all stabilizers and implements a non-trivial transfor- 
mation on the logical qubits, see [17] ■ Let us show that 
the subsystem toric code has distance d — L. Indeed, 
any error E as above must anti-commute with at least 
one of the logical operators X\, Zi,X 2 , Z 2 . Assume wlog 
that E anti-commutes with Z\. Let A' be any vertical 
line on the lattice (a horizontal translation of A) and 

z\ = YijeA' Zj- One can easily check that Z\Z\ coin- 
cides with the product of stabilizers Sp over all plaque- 
ttes p lying between A and A'. Since E commutes with 
all stabilizers, we conclude that E anti-commutes with 
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Z-y. But this means that E must act non-trivially on at 
least one qubit of A'. Since there are L non-overlapping 
choices of the line A', we conclude that E must have 
weight at least L. One can also easily check that trans- 
lating A by the half of the lattice period gives a logical 
operator of weight L equivalent to Z\ (modulo gauge op- 
erators). Hence the code has distance d = L. 

The first step in any error correction protocol based 
on stabilizer codes is the syndrome readout, that is, a 
non- destructive eigenvalue measurement of every stabi- 
lizer operator. To measure the 6-qubit stabilizers S p 
and Sp we shall take advantage of the gauge qubits 
and the identity Eq. The simplest syndrome read- 
out protocol consists of two steps: Step 1. Measure the 
eigenvalue of every X-type triangle operator G{T) and 
record the outcome A(T) = ±1. Step 2. Measure the 
eigenvalue of every Z-type triangle operator G(T) and 
record the outcome A(T) = ±1. Since any triangle op- 
erator commutes with stabilizers, the eigenvalue of any 
stabilizer remains unchanged throughout the execution 
of this protocol. Hence the eigenvalues of stabilizers 
S% and S* are given by A(S^) = \(T p SE )\(T p NW ) and 
A(S**) = \(T p sw )\{T p NE ), see Eq. g). In practice it may 
be advantageous to use 'interleaved^ protocols in which 
Steps 1,2 defined above are implemented in parallel, see 
Section [Vl for more details. 



III. TOPOLOGICAL QUANTUM ORDER 



Consider a Hamiltonian 



H 



(6) 



where the sum is over all triangles of the lattice. Re- 
call that G(T) are the 3-qubit triangle operators defined 
in Eq. 0. In this section we compute the entire eigen- 
value spectrum of H and show that on the torus H has 
a four-fold degenerate ground state separated from ex- 
cited states by a constant energy gap. Moreover, we shall 
construct a unitary locality preserving transformation U 
such that UHW can be regarded as the standard toric 
code Hamiltonian on the square lattice (with some irrele- 
vant ancillary qubits) . Thus the model defined in Eq. ([6| 
exhibits topological quantum order. 

Let us first compute eigenvalues of H. Since the sta- 
bilizers Sp , Sp commute with every term in H, we can 
assume that any eigenvector ip of H is also an eigenvector 
of any stabilizer, that is, ip = x p ip and Spip — z p i p for 
some syndromes x p ,z p = ±1. Using identities Eq. (2|3) 
one gets 

Hip = - ^(1 + x p )Xpi> + (1 + Zp)Z p tP 



where X pi Z p are the logical operators on the gauge qubit 
located at the plaquette p. Hence the restriction of H 



onto the sector with fixed syndromes x p , z p describes L 2 
non-interacting gauge qubits. 

Let eo(x p ,Zp) and ei(x p ,Zp) be the smallest and the 
largest eigenvalues of a gauge qubit p for a fixed syn- 
dromes x p , z p . A simple algebra shows that 
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To minimize the overall energy one has to choose x p = 
z p = 1 for all p. This shows that ground states of H 
belong to the trivial syndrome sector and the ground 
state energy is E a = —2\/2L 2 . The ground state is four- 
fold degenerate since the code has two logical qubits. 

Excitations of H fall into two categories. First, there 
are gauge excitations that are confined to the trivial syn- 
drome subspace x p — z p = 1. The energy cost of a single 
gauge excitation is A g = ei(l, 1) — eo(l, 1) = 4-\/2. A 
gauge excitation on a plaquette p can be created locally 
by a proper combination of operators X p and Z p . Hence 
gauge excitations do not carry any topological charge. 
Secondly, there are syndrome excitations that flip syn- 
drome bits x p and z p . The energy cost of a single syn- 
drome excitation is A s = e (l, —1) — eo(l, 1) = 2(v2 — 1). 
It corresponds to flipping x p (or z p ) on any plaquette p. 
A single syndrome excitation however cannot be created 
locally due to the constraints J7 p x p = lip Z P — 1, see the 
previous section. It means that syndrome excitations can 
only be created in pairs. Each pair costs energy 2A S . 

We can now show that the Hamiltonian of Eq. ^ is lo- 
cally equivalent to Kitaev's toric code. Consider a quan- 
tum circuit U shown on Fig. [3] It consists of four rounds 
of CNOT gates, U = U^U^U^U^, where E/« is a 
tensor product of L 2 disjoint CNOT gates labeled by j on 
Fig. [3] Note that U is a locality preserving transforma- 
tion, that is, the Heisenberg evolution of any observable 
O — > UOTP can only enlarge the support of O by a few 
units of length. Such transformations do not change any 
topological features of the model [53]. A simple alge- 
bra shows that the transformed stabilizers USpW = A p 
and USpW = B p coincide with the star and plaquette 
operators of the Kitaev's toric code on a tilted square 
lattice, see Fig. [3] Furthermore, the transformed gauge 
generators UX p W = J* and U Z p W = Jp become one- 
qubit Pauli operators X and Z respectively acting on the 
qubit located at the bottom edge of p, see Fig. [3| Using 
identities Eqs. (|2|3|) we arrive at 



H' = UHU^ 



p 



Jp Up. 



The same arguments as above show that ground states 
of H' are defined by equations A p ipa — B p ipo — ipo and 
(Jp + Jp)ipa = v^V'o f° r au plaquettes p. Thus any 
ground state of H' must have a form ipQ = ipt op <8> ipanc. 
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where iptop is a ground state of the Kitaev's toric code 
on the tilted square lattice while 4>anc is a tensor product 
of one-qubit ancillary states located on horizontal edges 
of the original lattice. Such ancillary unentangled states 
clearly have no effect on topological features of the model. 
We conclude that that the Hamiltonian Eq. ^ is in the 
same topological phase as the toric code model. The 
exact solvability of the model clearly extends to a more 
general Hamiltonian H — — X)t 9tG(T), where gx are 
arbitrary coefficients. 

Consider now a modified Hamiltonian 

H" = -J2 Jp + Jp +A P + B p . 

p 

Note that H' and H" have the same ground subspace 
and H" coincides with the ordinary toric code Hamilto- 
nian p] if one ignores the ancillary qubits. For any state 
<f> orthogonal to the common ground subspace of H' , H" 
and for any parameter < t < 1 one has 

(0|(1 - t)H' + tH"\<f>) > (1 - t)A' + tA" 

where A' = 4(\/2 — 1) and A" = 2 are the energy gaps of 
H' and H" respectively. It follows that the Hamiltonian 
(1 - t)H' + tH" has energy gap at least min (A', A") = 
4(\/2 - 1) for all < t < 1. Hence we can continuously 
deform H' to H" without closing the gap. 

The decoupling circuit of Fig. [3] illustrates a feature of 
a more general family of 2D translationally-invariant sub- 
system codes described in [50]: they can be mapped to 
(one or more copies of) the standard toric code with ad- 
ditional decoupled qubits by a constant-depth quantum 
circuit. 




FIG. 3: Decoupling circuit. Applying four rounds of CNOT 
gates as shown on the left in a translational invariant way 
transforms the gauge operators and stabilizers as shown on 
the right. The stabilizer generators become those of Kitaev's 
toric code on a tilted square lattice, while the extra gauge 
generators are mapped to single-qubit Pauli operators acting 
on ancillary qubits (horizontal edges). Thus, the circuit has 
the effect of decoupling the gauge operators from the toric 
code. 

We note that in general the Hamiltonian obtained by 
the sum of the gauge operators of a topological subsystem 



code does not necessarily produce topological order. The 
peculiar feature of the present model is that the canonical 
logical Pauli operators on the gauge qubits X p and Z p are 
local. In contrast, it was shown in |25 that the subsystem 
color code can be obtained from multiple copies of the 
ordinary color code by gauging out both local and non- 
local logical operators. In particular, some of the logical 
operators that are gauged out carry topological charge. 
The present analysis does not apply to such models. 

IV. SUBSYSTEM SURFACE CODES 

We can now describe a subsystem version of the sim- 
plest surface code on a planar lattice with two rough and 
smooth boundaries that encodes one logical qubit [2]. 
Now the lattice has open boundary conditions. A lat- 
tice of size L x L has (L + l) 2 vertices, 2L(L + 1) edges, 
and L 2 plaquettes. As before, code qubits are placed at 
vertices and centers of edges, so the total number of code 
qubits is n = (L + l) 2 + 2L(L + 1) = 3L 2 +4L + 1. 

For every edge e lying on the boundary of the lattice 
define a weight-2 Pauli operator S e of type XX or ZZ 
as shown on Fig. |4j Let S be the group generated by 
wcight-6 operators S p , S p defined earlier and all wcight- 
2 operators S e associated with boundary edges. One can 
easily check that S is an abelian group with s = 2L 2 +4L 
independent generators. Therefore it defines a stabilizer 
code encoding k' = n — s = L 2 + I qubits. 



z z 




XX z z 



FIG. 4: Subsystem version of the surface code. 

Since the additional weight-2 stabilizers S e lying on the 
boundary commute with all triangle operators G(T), we 
can use Eq. ^ to define g = L 2 gauge qubits associated 
with plaquettes of the lattice. Logical operators on the 
remaining k = k' — g = 1 logical qubit can be chosen as 
X\ and Z\, see Eq. Q, that is, a horizontal line of X'b 
and a vertical line of Z's. Note that a vertical line of X's 
and a horizontal line of Z's are no longer logical operators 



() 



because they anti-commute with some of the boundary 
stabilizers S e . The same arguments as above imply that 
the code has minimum distance d — L. An extension to 
a planar lattice with punctured holes encoding multiple 
logical qubits is sketched on Fig. [5] 

z z z z 




z z z z 



FIG. 5: Subsystem version of the surface code with two punc- 
tured holes. Each pair of holes encodes one logical qubit. The 
external boundary of the lattice generates the third logical 
qubit. Thick red and blue lines show the logical operators. 



V. ERROR CORRECTION WITH NOISELESS 
SYNDROMES 

In this section we propose an error correction proto- 
col for the idealized setting where the syndrome readout 
is noiseless. For simplicity, we shall first focus on the 
subsystem toric code defined in Section [TT| It suffices to 
construct a protocol for correcting errors of A-type (bit 
flips). Due to the code symmetry, the same protocol can 
be applied to Z-type errors. Let E be an unknown Pauli 
error of X-type. We will say that E creates a defect 
at a plaquette p iff E anti-commutcs with the stabilizer 
S? . The syndrome measurement reveals a configuration 
of defects created by E. The key observation is that any 
single-qubit A-error creates exactly two defects. Indeed, 
an A error on a vertical or horizontal edge e creates de- 
fects at the two plaquettes adjacent to e. An X error at 
any vertex u creates defects at the two plaquettes lying 
in the north-west and the south-east quadrants of u. The 
relationship between single-qubit A-errors and the corre- 
sponding pairs of defects can be captured by introducing 
a virtual lattice A that consists of virtual vertices and 
virtual edges. A virtual vertex p represents a stabilizer 
Sp (plaquette p of the original lattice), while a virtual 
edge represents a pair of defects that can be created by a 
single-qubit X error. One can easily check that A is the 
regular triangular lattice, see Fig. [6j 

Furthermore, for each virtual edge e = (p, q) there is 
only one A-error creating a pair of defects at p and q. For 



any virtual vertex p let 6(j>) be the set of virtual edges 
incident to p. Then 

e€<5(p) 

that is, Z-type stabilizers can be regarded as star opera- 
tors of the standard toric code defined on the triangular 
lattice. Note that a closed loop on the virtual lattice en- 
closing any triangular face is an A-type triangle operator 
G(T), while non-contractible closed loops correspond to 
logical operators. 




FIG. 6: The virtual lattice A (solid lines) describing correc- 
tion of X-type errors for the subsystem toric code of Fig. [T] 
Opposite sides of the lattice must be identified. The original 
square lattice is shown by dashed lines. Each virtual edge 
represents a code qubit. Stabilizers of Z-type Sp correspond 
to stars on the virtual lattice (solid blue lines). Triangle oper- 
ators of X-type correspond to triangular faces of the virtual 
lattice (not shown). Error correction amounts to finding the 
minimum weight matching of defects on the virtual lattice. 
The virtual lattice describing correction of Z-type errors is 
obtained from A by the 90° rotation. 

Assuming that errors on different qubits are indepen- 
dent and have the same rate p, the most likely error E* 
consistent with the observed syndrome coincides with the 
minimum weight matching of defects on the virtual lat- 
tice. The latter can be found efficiently using the Ed- 
monds's algorithm, see [3j for details. Choosing E* as a 
correction operator always returns the system back to the 
codespace C. The overall evolution of the system is de- 
scribed by an operator EE* which has trivial syndrome 
and thus can be viewed as a linear combination of closed 
loops on the virtual lattice with Z2 coefficients. Error 
correction is successful iff EE* acts non-trivially only on 
the gauge qubits, that is, EE* is a product of A-type tri- 
angle operators G(T). Equivalently, EE* must represent 
the trivial cycle in the homology group i?i(A, Z2). 

As was argued in [3l [HJ [14] , the optimal error correc- 
tion strategy amounts to finding the most likely equiv- 
alence class of errors consistent with the observed syn- 
drome rather than the most likely error. More specifi- 
cally, let G be the gauge group generated by the triangle 
operators G(T) of A-type. Since logical operators act- 
ing on gauge qubits are irrelevant, all errors in the coset 
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E ■ Q must be considered equivalent. Note that there are 
only four cosets of Q consistent with the observed syn- 
drome, namely, E ■ G, EX 1 ■ Q, EX 2 ■ Q, and EX 1 X 2 • Q. 
As was shown in [5J [T3], the probability of each coset 
can be expressed as the partition function of the random 
bond ± J Ising model. In our case Ising spins reside on 
triangular faces of the virtual lattice, anti-ferromagnetic 
bonds correspond to virtual edges that belong to E, and 
the inverse temperature (3 is determined by the Nishimori 
condition e~ 2/3,7 = j^-. The threshold error rate p c for 
the optimal decoding coincides with the density of anti- 
ferromagnetic bonds at the phase transition point [3]. 
The latter has been recently computed by Queiroz 
who found p c m 7%. The analogous threshold error rate 
for the standard toric code is known to be approximately 
11%, see [3]. 

One can similarly construct the virtual lattice for the 
subsystem surface code, see Fig. [7] The only difference 
is that now defects can be matched either to each other, 
or to one of the boundaries. 




FIG. 7: The virtual lattice A (solid lines) describing correc- 
tion of X-type errors for the subsystem surface code of Fig. [4] 
Diamonds and open circles represent stabilizers Z-type stabi- 
lizers Sp and S e respectively. The virtual lattice describing 
correction of Z-type errors is obtained from A by the 90° 
rotation. 



VI. ERROR CORRECTION FOR THE 
CIRCUIT-BASED ERROR MODEL 

Let us now consider more realistic settings when the 
syndrome information itself may contain errors. We as- 
sume that the library of elementary operations supported 
by the quantum hardware includes CNOT gates between 
nearest-neighbor qubits, single-qubit measurements in 
X- or Z-basis, and preparation of single-qubit ancillary 
states |0) and |+). Our error correction protocol will be 
defined as a sequence of rounds, where at each round any 
qubit can participate in one elementary operation. Wc 
assume that each elementary operation is noisy, so it can 
fail with a probability p that we call an error rate. More 
precisely, our error model, borrowed from [5], is defined 



as follows. 

• A noisy X or Z measurement is the ideal measure- 
ment in which the outcome is flipped with proba- 
bility p. 

• A noisy |0) or |+) ancilla preparation returns the 
correct state with probability 1—p and the orthog- 
onal state with probability p. 

• A noisy CNOT gate is the ideal CNOT followed by 
an error (1 — p)Id + pD, where Id is the identity 
map and T> is the fully depolarizing two-qubit map 
applying one of 16 two-qubit Pauli operators with 
probability 1/16 each. 

We assume that once a qubit has been measured, its state 
is unknown. To use such a qubit again, it must be explic- 
itly initialized using the noisy preparation defined above. 
We do not need to dehne memory errors because no qubit 
will be idle at any round of our protocol. 

In order to measure eigenvalue of individual triangle 
operators ZiZjZk and X p X q X r we shall use quantum 
circuits shown on Fig. [8| Measuring a single triangle 
operator requires one ancillary qubit and five rounds. 
Similar circuits with one extra CNOT gate were used 
in fault-tolerant protocols based on the standard surface 
code [3J IB] , where one has to measure four-qubit plaque- 
tte and star operators Z® 4 and X® 4 . 
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FIG. 8: Quantum circuits for measuring the eigenvalue of tri- 
angle operators ZiZjZu (top) and X p X q X r (bottom). The 
circuits use one ancillary qubit. A single Z (X) error on the 
ancilla a z (a x ) can propagate via CNOTs to at most one Z 
(X) error on code qubits modulo gauge operators. A single 
X (Z) error on the ancilla a 2 (a x ) results in a faulty measure- 
ment outcome. 

We begin by highlighting strengths and weaknesses of 
the subsystem and the standard surface code. The key 
advantage of the SSC is a limited propagation of errors 
by the CNOT gates. Consider, for example, a single Z er- 
ror on the ancillary qubit a z in the circuit measuring the 
triangle operator ZiZjZk, see Fig. [8] (top). Depending 
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on the round at which the error has occurred, it propa- 
gates to one of the errors Zj, Z^Zj, ZiZjZk on the code 
qubits. Multiplying the last two errors by the triangle 
operator Z^ZjZfc leaves only singlc-qubit errors Zj, Zk, 
and the identity error. It shows that a single Z error on 
the ancilla a z can lead to at most one Z error on the 
code qubits, modulo gauge operators. An X error on a z 
cannot propagate to code qubits, so its only effect is flip- 
ping the measurement outcome which results in a faulty 
syndrome bit. Since each stabilizer is composed of two 
disjoint triangle operators, see Fig. |4j each syndrome bit 
effectively accumulates errors from ten rounds. For com- 
parison, the standard surface code requires six rounds to 
measure a single syndrome bit, however a single error on 
the ancilla can feed back to two errors on code qubits 
(such double errors were referred to as 'horizontal hooks' 
in [3]). This shows that neither of the two codes offers 
an obvious advantage compared with the other. 

Let us now discuss our syndrome readout circuit in 
more details. Since individual syndromes can no longer 
be trusted, we shall repeat syndrome measurements T 
times for some T> 1. We choose T = L in all numeri- 
cal simulations. Error correction is deemed successful if 
the accumulated error E on the code qubits that results 
from T noisy syndrome measurements can be corrected 
based on the full observed syndrome history and one final 
syndrome readout which we assume to be noiseless. (In 
practice the final readout involves measuring each code 
qubit in |0) or |+) basis. Outcomes of such measurement 
determine the syndrome of Z-type or A-type stabilizers 
respectively. We can assume that single-qubit measure- 
ments are noiseless by absorbing measurement errors into 
memory errors that occurred one round earlier.) 

Repeating the circuits shown on Fig. [8] cyclically L 
times would naively require 5L rounds. We can reduce 
the required number of rounds to 4L by introducing two 
ancillary qubits for each triangle. One of them serves as 
the ancilla a z or a x shown on Fig. [8] The purpose of the 
second ancilla is to enable offline preparation of |0) or |+) 
states which can be performed in the same round as the 
measurement of the main ancilla. To simplify notations, 
we will only show one ancilla per triangle and assume 
that this ancilla is initialized in the |0) or |+) state at 
the end of each measurement round (with a duly added 
noise) . 

The readout circuit will be chosen such that any fixed 
triangle alternates between three gate rounds and one 
measurement round in a cyclic fashion. It can be repre- 
sented by a local readout schedule 

• • • MGiG 2 G 3 MGiG 2 G 3 M ■ ■ ■ (7) 

where M is either X-type or Z-type measurement on 
the ancilla while G*i , G 2 , G 3 are CNOT gates coupling 
the ancilla and the code qubits (the ancilla is control for 
X-type triangles and target for Z-type triangles). Time 
flows from the left to the right. Note that each trian- 
gle has 24 different choices of its local schedule. Indeed, 
there are 6 choices of the order in which the ancilla is 



coupled to the code qubits and 4 choices of the round at 
which the triangle is measured. Local schedules chosen 
at different triangles must be consistent with each other, 
such that at any round any qubit participates in at most 
one operation. 

We shall focus on schedules which are periodic both 
in space and time. Hence the entire readout circuit is 
completely specified by local schedules inside a 3D el- 
ementary cell which consists of four rounds labeled as 
0, 1, 2, 3 and four triangles of type NE, SE, SW, NW lo- 
cated at some fixed plaquctte p. We begin by observing 
that a consistent schedule cannot have a round at which 
every triangle applies a CNOT. Indeed, this would define 
a matching between code qubits and triangles. However, 
the lattice has 4L 2 triangles and only 3L 2 code qubits. 
This observation shows that at every round all triangles 
of some type have to be measured. Furthermore, it is 
natural to demand that if some pair of triangles form a 
stabilizer, these triangles must be measured in two con- 
secutive rounds (otherwise the corresponding syndrome 
bit would accumulate too much errors). In other words, 
we would like to measure all triangles of X-type in two 
consecutive rounds and all triangles of Z-type in the two 
remaining consecutive rounds. Our choice of the mea- 
surement rounds satisfying these conditions is shown on 
Fig. [9] Note that all X-type triangles are measured at 
rounds 3, 0, while all Z-type triangles are measured at 
rounds 1,2. 




Measurement rounds 

FIG. 9: The numbers indicate rounds (modulo four) at which 
the ancillas assigned to each triangle have to be measured. 

It remains to schedule CNOT gates. We will say that 
a schedule is correct iff for each triangle one can move 
all gates forward in time towards the next measurement. 
Here moving a gate is allowed as long as it commutes 
with other gates. A correct schedule faithfully simulates 
the simple syndrome extraction routine described in Sec- 
tion [H] since after moving all gates towards the next mea- 
surement all X-type stabilizers are measured at rounds 
3, 0, while all Z-type stabilizers are measured at rounds 
1, 2. We shall look for a schedule which is correct and in- 
variant under the exchange of X and Z (modulo lattice 
symmetries and time translations). The latter require- 
ment is rather natural since our error model does not 
have a bias towards X or Z errors. 
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To derive sufficient conditions for correctness let us 
introduce some terminology. Consider any triangle T and 
its local schedule, see Eq. ([7| . We shall refer to the gates 
G\ and G3 as the first gate and the last gate for the chosen 
triangle. If some pair of triangles T and T' are measured 
at two subsequent rounds j and j + 1 respectively, we 
will say that T is the leading triangle and T' is the tailing 
triangle. If T and T' are measured two rounds apart, the 
choice of the leading and the tailing triangle is arbitrary. 

Lemma 1. A schedule of CNOTs is correct if for any X- 
type and any Z-type triangle at least one of the following 
is true: 

• The two triangles are disjoint, 

• The two triangles are measured two rounds apart, 

• The last gate of the leading triangle commutes with 
the first gate of the tailing triangle. 

Proof. Suppose T x and T z are X-type and Z-type trian- 
gles respectively. If T x and T z are measured two rounds 
apart, their combined local schedules can be represented 
by a diagram 





M x 


Gf 


G% 


Gf 


M x 


Gf 


Gf 






GI 


Gf 


M z 


Gf 


G| 


Gf 


M z 





The gates Gf and Gf must be disjoint. Similarly, the 
gates Gf and Gf must be disjoint. Hence we can de- 
form the diagram by moving , one round forward 
and moving Gf, Gf one round backward obtaining an 
equivalent circuit: 





M x 


/~ix r~<x /~*x 
( - r l ( - r 2 ( - r 3 


M x 


t~ix r~<x r~<x 
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M z 
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We can further deform the circuit by moving each mea- 
surement backwards towards the next gate. 

Suppose now that T x and T z are measured in two sub- 
sequent rounds such that T x is the leading and T z is 
the tailing (the opposite case is completely analogous). 
Then their combined schedules can be represented by a 
diagram 





M x 


Gf 


Gf 


Gf 


M x 


Gf 






G\ 


M z 


Gf 


G z 2 


G\ 


M z 





By assumption, the gates Gf and Gf commute. Hence we 
can deform the diagram by moving Gf one round back- 
ward and moving Gf one round forward. In addition, we 
can move Gf one round forward and move Gf one round 
backward. It gives an equivalent circuit: 





M x 




(~ix /~*x r~<x 
( - r l ( - r 2 ( - T 3 




M x 




r~ix r^x r^x 








M z 




r^z t~iz /~*z 
( - T l ( - r 2 l - r 3 




M z 







We can further deform the circuit by moving each mea- 
surement backwards towards the next gate. 

After the deformation each triangle applies the en- 
tire sequence G1G2G3M in a single round which is two 
rounds apart from the measurement round in the origi- 
nal schedule. We can now move the sequence G1G2G3M 
two rounds forward for each triangle simultaneously un- 
til X-type triangles apply Gf Gf Gf M at rounds 0, 3 and 
Z-type triangles apply GfGfGfM at rounds 1,2. This 
shows that the original schedule is correct. □ 




Gate rounds 



FIG. 10: Example of a correct syndrome extraction sched- 
ule with four rounds labeled by 0,1,2,3 repeated cyclically. 
The numbers assigned to vertices of each triangle indicate the 
rounds at which the code qubits comprising a triangle are cou- 
pled to the ancilla by CNOT gates. Measurement rounds are 
shown on Fig. [9] One can check that any pair of a>type and 
z-type triangles satisfies conditions of Lemma [I] 

Our choice of a CNOT schedule is shown on Fig. [10] 
One can easily check that it satisfies conditions of 
Lemma [T] It remains to define the classical post- 
processing step that extracts the syndromes from the 
measured eigenvalues of triangle operators. For any inte- 
ger t G [0, L — 1] and a plaquette p we define a syndrome 
bit Sp (t) as a product of eigenvalues of Z-type triangles 
SE and NW located at the plaquette p that were mea- 
sured at rounds At + 1 and At + 2. Similarly, we define a 
syndrome bit Sp (t) as a product of eigenvalues of X-type 
triangles SW and NE located at the plaquette p that were 
measured at rounds At + 3 and At + A. Hence each syn- 
drome bit combines eigenvalues of two triangle operators 
measured in two consecutive rounds. 

Let us now move to the error correction protocol that 
takes as input the syndrome information and outputs a 
correcting Pauli operator E* acting on the code qubits. It 
mostly follows 3, SJ [2T] . Our protocol deals with X-type 
and Z-type errors independently. It should be noted that 
the schedule shown on Figs. |9|10| is invariant under the 
horizontal reflection of the lattice and shifting the time 
by two rounds. Since the horizontal reflection exchanges 
X-type and Z-type triangles, it suffices to analyze X- 
type errors. 

Let us introduce a 3D virtual lattice A that consists of 
virtual vertices and virtual edges. A virtual vertex is a 
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pair (p, t) , where p is a plaquette of the 2D code lattice 
and t £ [0, L— 1] is the discrete time. We shall say that a 
virtual vertex u = {p, t) has a defect iff the syndrome bits 
Sp (t) and Sp (i + 1) are different. Hence the full syndrome 
history can be regarded as a configuration of defects on 
the virtual lattice. 

We begin by considering configurations of defects cre- 
ated by a single fault in the readout circuit. Here a single 
fault includes one of the following possibilities: 

• Wrong measurement outcome on some ancilla, 

• Wrong ancilla preparation, 

• One of the errors IX, XI, or XX inserted after 
some CNOT gate. 

In other words, a single fault is any event that can occur 
with probability Q(p) in the limit p — > (recall that we 
only keep track of A- type errors). In order to define 
virtual edges we need the following observation. 

Lemma 2. Any single fault in the readout circuit creates 
either or 2 defects on the virtual lattice. 

Proof. A measurement error on the ancilla a z at a pla- 
quette p creates two defects at virtual vertices (p, t) and 
(p,t + 1) for some t. Measurement errors on ancillas a x 
create no defects since we ignore A-type syndromes. An- 
cilla preparation error can be regarded as an X-type error 
for ancillas a z and Z-type error for ancillas of . Such an 
error can be propagated forward in time without feeding 
back to the code qubits because a z is always a target 
qubit and a x is always a control qubit for any CNOT 
gate, see Fig. [8] Hence ancilla preparation errors are 
equivalent to the measurement errors. By the same rea- 
son, an A-error on a z caused by any CNOT gate is equiv- 
alent to a measurement error. 

Consider now a single-qubit A error on some code 
qubit. If the syndrome were measured on all plaquettes 
directly after the error, one would observe non-trivial 
syndromes s^ , s^ at some pair of plaquettes p,q, see 
Section [U Since there are no other errors in the readout 
circuit, it will faithfully simulate the ideal syndrome mea- 
surements, that is, the syndrome s^ (t) will change from 1 
to —1 for some step t p and the syndrome s q (t) will change 
from 1 to —1 for some step t q . This produces a pair of 
defects at virtual sites (p,t p ) and (q,t q ). (More detailed 
analysis shows that either t p — t q or t p — t q ± 1). An er- 
ror AA that occurred after a CNOT gate is equivalent to 
a single A error on the control qubit that occurred before 
this CNOT. The only remaining case is a single A error 
on the ancilla of . It can be propagated forward or back- 
ward towards the nearest admeasurement. Such prop- 
agation feeds back at most one A error to code qubits. 
This is the case that we have already explored. □ 

A trivial corollary of the lemma is that the number of 
defects on the virtual lattice is always even. 

We connect a pair of virtual vertices u, v by a virtual 
edge iff there a single fault in the readout circuit capable 



of creating a pair of defects at u and v. A more detailed 
analysis shows that the virtual lattice has seven types of 
edges (not counting the orientation). The table below 
shows all 14 neighbors v of some fixed virtual vertex u = 
(x, y, t) as well as the total number of single faults of each 
type in the readout circuit that create defects at u and v. 
For brevity we refer to single faults IX, XI, XX as G- 
faults (gate faults), while measurement and preparation 
single faults are referred to as M-faults and P-faults. 



Neighbor of (x,y,t) 


G-faults 


M-faults 


P-faults 


Prior 


(x,y,t± 1) 


6 


2 


2 


Up/2 


(x± l,y,t) 


8 








2p 


(i±l,l/Tl,t) 


8 








2p 


[x,y± l,t) 


4 








P 


(x,y±l,t±l) 


2 








p/2 


(xt l,j/,t±l) 


2 








p/2 


(iTl,l/±l,t±l) 


2 








p/2 



Note that the space-like virtual edges (those for which u 
and v have the same t coordinate) correspond to edges of 
the 2D virtual lattice defined in Section [V] If one ignores 
the t coordinate, space-like virtual edges correspond to 
the code qubits. In particular, pairs of defects located on 
space-like virtual edges can be viewed as memory errors. 
On the other hand, pairs of defects located on time-like 
virtual edges(those for which u and v have the same x,y 
coordinates) can be viewed as syndrome measurement 
errors. The remaining virtual edges represent various 
combinations of memory errors and measurement errors. 

For every virtual edge e we define a prior p e as the 
probability to observe a pair of defects at the endpoints 
of e. Taking into account that any single G-fault has 
probability p/A, while a single M-fault and a single P- 
fault have probability p, one arrives at the priors listed 
in the table. 

We shall choose the correction operator E* by pre- 
tending that the creation of defect pairs on different vir- 
tual edges are independent events. Then the most likely 
combination of memory errors and measurement errors 
consistent with the observed configuration of defects co- 
incides with the minimum weight matching of defects on 
the virtual lattice, where an edge e is assigned a weight 
w e ~ log (l/p e ). The minimum weight matching M can 
be found efficiently using the Edmonds's algorithm. Fi- 
nally, we choose the correction operator E* as the prod- 
uct of all memory errors that appear in M. In order to 
decide whether the error correction is successful we com- 
pare E* with the accumulated error E on the code qubits 
generated by the syndrome readout circuit. The results 
of our Monte Carlo simulation are shown on Figs. |11|12| 
It indicates that the threshold error rate for the circuit- 
based error model is p c ~ 0.6%. 
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FIG. 11: Probability of the error correction failure for the sub- 
system toric code under the circuit-based error model. For a 
lattice of size L the syndrome measurement has been repeated 
L times. Each data point was obtained using 10 4 - 10 6 Monte 
Carlo trials. The simulation was performed only for X-type 
errors (Z-type errors on the reflected lattice). 
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FIG. 12: Scaling of the error correction failure probability for 
different error rates below the threshold. Error bars represent 
the statistical error. Achieving the failure probability 10 -6 at 
the error rate p = 0.5% would require L « 100. 



VII. DIRECT 3-QUBIT PARITY 
MEASUREMENTS 

In the previous section we estimated the threshold in 
a model where the triangle operators are measured by 
performing CNOT gates between the code qubits and an 
ancillary qubit. In this section, we will study a differ- 
ent model where measurements are performed directly. 
For some physical devices, it is possible to directly probe 
a multi-qubit parity operator without the need to use 
an ancillary qubit — the probe itself is used as a medi- 
ator that stores and accumulates the multi-qubit corre- 
lations. Direct two-qubit measurements have been re- 
alized in circuit quantum electrodynamics [28-31], and 



there are proposals [33] to turn these two-qubit measure- 
ments into parity measurements. These direct measure- 
ment schemes require a separate threshold analysis be- 
cause they have different noise models and propagation. 

We will focus on the recent proposal of DiVincenzo 
and Solgun |21j that realizes a three-qubit parity mea- 
surement ZZZ by capacitively coupling three Josephson- 
junction qubits to two transmission-line resonators. In 
the dispersive regime — where the difference between the 
resonant frequency of the transmission lines and the qubit 
transition frequency is much larger than the coupling 
strength — the transmission-line will pick-up a qubit- 
state-dependcnt frequency shift (Stark shift). When a 
near-resonant frequency probe signal is sent through one 
of the transmission lines, it picks up a phase that depends 
on the resonant frequency. Thus, the state-dependent 
resonant frequency shift will imprint a phase shift on 
the probe pulse that depends on the state of all three 
qubits. With an appropriate choice of parameters (qubit- 
transmission line detuning, coupling strength, probe sig- 
nal frequency), the probe signals sent through the two 
transmission lines can be measured interferometrically 
to reveal information only about the parity of the three 
qubits, all other information about the qubit state im- 
printed on the probe signals being erased by the interfer- 
ometric measurement. One important advantage of this 
measurement scheme is that a single qubit can partici- 
pate to two distinct parity measurements simultaneously 
with an appropriate arrangement of transmission lines. 

The main source of errors in this measurement is de- 
phasing caused by the finite bandwidth of the probe pulse 
|21) . In an ideal parity measurement, the coherence be- 
tween two computational basis states of the three-qubit 
systems \x) and \y) would be totally suppressed when x 
and y have different parities, and unaffected otherwise. 
The finite bandwidth of the probe pulse will cause de- 
phasing between states of a given parity, and incomplete 
dephasing between states of distinct parity. Moreover, 
these errors are otherwise uniform, they do not depend, 
e.g., on the Hamming distance \x— y\ between the compu- 
tational basis states. This is characteristic of a collective 
noise model, where multi-qubit errors are as likely as sin- 
gle qubit errors. In contrast, when qubits are subject to 
independent noise, dephasing would increase with Ham- 
ming distance between the states. 

Parity measurements in the conjugate basis are re- 
quired to measure the X-type triangle operators. These 
can be realized by rotating the qubits prior to sending the 
probe signal in the transmission line. Single qubit rota- 
tions are very fast and accurate in this architecture |33j . 
Nonetheless, they propagate errors, and can interchange 
X-type and Z-type errors. Based on these considera- 
tions, we will model noisy measurement of X-type and 
Z-type triangles in the following way: 

• A noisy XXX or ZZZ measurement is modeled by 
a perfect even/odd subspace projection, followed by 
an error (l—p)Id+pD, where Id is the identity map 
and T> is the fully depolarizing three-qubit map ap- 
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plying one of 64 three-qubit Pauli operators with 
probability 1/64 each. 

• The measurement outcome is flipped with proba- 
bility p. 

With this model, the syndrome extraction cycle requires 
only two rounds: one to measure all X-type triangles 
and one to measure all Z-type triangles. This implies 
that some qubits participate to two simultaneous mea- 
surements. As mentioned above, this is not a problem 
physically so far as a single qubit can be coupled to mul- 
tiple transmission lines, which has already been demon- 
strated experimentally [34 . Moreover, two noisy parity 
measurements in the same basis always commute with 
our choice of noise model. 

We have simulated fault-tolerant error correction of 
the subsystem toric code using direct parity measure- 
ment with the noise model described above, our results 
are presented on Fig. |13| Since errors are correlated 
in this model, we have opted for the renormalization 
group (RG) decoding algorithm proposed in [HI [12] . In- 
deed, Edmonds's minimum weight matching algorithm 
assumes an independent noise model and consequently 
yields a lower threshold in the presence of correlations. 
The RG decoder can exactly account for some of these 
correlations (those that are contained within a RG unit 
cell). Additional correlations can be approximated by 
updating the error prior on each RG unit cell using a 
belief propagation decoder [35] . The RG decoding is ex- 
ecuted using these updated error priors. Following the 
proposal of [5D] , we map the code and its (updated) noise 
model onto the standard toric code on which we execute 
RG [35] ■ While Refs. PUCES] assumed noiseless syn- 
drome measurements, the decoding algorithm can easily 
be extended to noisy syndrome by renormalizing a 3D 
lattice l36l. 
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FIG. 13: Decoding error probability as a function of the phys- 
ical error rate p for various lattice sizes L. The lines are ob- 
tained using the fitting functions a + b(p — pth)L x ' u + c[(p — 
Pth)L 1 / v ] 2 , which produces a threshold pth ~ 0.968% and a 
critical exponent v w 1.36. 



The results shown on Fig. 13 indicate a threshold value 
of roughly 0.97%. This value should be seen as lower 
bound to the true threshold of this code, which may well 
be above 1%. Indeed, the results presented here were 
obtained from a unit cell of dimension 2x2x1: two 
of the three space-time dimensions are renormalized at 
each iteration, and by rotating the unit cell at each iter- 
ation we obtain a renormalization of the entire space- 
time lattice by a factor 4 after 3 RG rounds. Based 
on our experiments, larger RG unit cells produce higher 
thresholds because they make use of more correlations 
existing in the noise model. The decoding complexity 
scales exponentially with the unit cell size however, so 
in practice we are limited to relatively small cells. Fur- 
thermore, while the simulated syndrome extraction pro- 
tocol involved measurements of both X-type and Z-type 
triangles, error correction has been performed only for 
X-type errors. Applying the same error correction algo- 
rithm independently to X-type and Z-type errors would 
result in the same error threshold [35]. We expect that 
more sophisticated decoders taking into account correla- 
tions between X and Z errors could achieve higher error 
thresholds. 

Finally, we note that direct parity measurements of 
weight-four operators such as ZZZZ and XXXX can 
be realized similarly [3T]. This could be used to imple- 
ment Kitaev's toric code. However, it is not reasonable 
to assume that the noise rate p is independent of the 
weight w of the operator being measured. Thus, one 
needs to work out this dependence p(w) from physical 
considerations before comparing thresholds of different 
codes obtained from direct parity measurements. 



VIII. CONCLUSION 

We have presented a subsystem version of Kitaev's sur- 
face code. The main features of our code is that it re- 
quires only 3-qubit parity measurements and its stabilizer 
generators have weight 6. Minimizing the weight of the 
parity measurements is helpful as it simplifies the mea- 
surement procedure, while minimizing the weight of the 
stabilizer generators is also desirable since it makes syn- 
dromes more reliable. In contrast to our code, the stan- 
dard toric code requires weight 4 parity measurements 
and has weight 4 stabilizer generators. The subsystem 
color codes require only weight 2 parity measurements, 
but have stabilizer generators of weight up to 18. Thus, 
based only on these considerations, it is not clear how 
the threshold of these various codes should compare. 

Our numerics show that in the circuit based model, our 
code has a threshold (0.6%) which is almost an order of 
magnitude larger than the one of the color code (0.08%) 
[37], and a bit more than half that of the standard toric 
code (0.9%) [5]. Motivated by the recent work of DiVin- 
cenzo and Solgun [21], we have also considered a setting 
where parity measurements can be implemented directly 
and found a threshold of 0.97%. This value cannot be 
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compared directly to the thresholds reported above since 
it is based on a substantially different noise model, and 
additional physical considerations must be taken into ac- 
count before comparing. 

We have shown that the new subsystem toric code gives 
rise to an exactly solvable spin Hamiltonian with 3-qubit 
interactions and topologically ordered ground state which 
is locally equivalent to the standard toric code. 
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